On the validity of the method of reduction of dimensionality: area of contact, average 

interfacial separation and contact stiffness 
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It has recently been suggested that many contact mechanics problems between solids can be 
accurately studied by mapping the problem on an effective one dimensional (ID) elastic foundation 
model. Using this ID mapping we calculate the contact area and the average interfacial separation 
between elastic solids with nominally flat but randomly rough surfaces. We show, by comparison to 
exact numerical results, that the ID mapping method fails even qualitatively. We also calculate the 
normal interfacial stiffness K and compare it with the result of an analytical study. We attribute 
the failure of the elastic foundation model to the neglect of the long-range elastic coupling between 
the asperity contact regions. 



1 Introduction 

The calculation of the stress and displacement field 
resulting from the contact between elastic solids with 
rough surfaces is a very complex problem, in part due 
to the many length scales usually involved, and also 
because of the long-range elastic coupling between the 
contact regions. For this reason simplifying approaches 
are very important. However, most analytical theories, 
such as the Greenwood- Williamson contact mechanics 
theoryjlj], or theories based on the elastic foundation 
modelQ (see Fig. Q]), neglect the elastic coupling be- 
tween asperity contact regions. It has recently been 
shown by exact numerical studies that the neglect of 
the elastic coupling results in qualitatively wrong con- 
tact topography 3f, and even the relation between the 
contact force and the area of contact is incorrectly de- 
scribed using this approach Q. In Ref. [H, [(| it was also 
shown that the contact stress-stress correlation function 
(in wavevector space) (<7(q)(7(-q)) ~ q~ a , where in the 
overlap model a = 2 + H (where H is the Hurst exponent), 
while including the long-range elastic coupling a = 1 + H . 

In a series of papers, Popov and coworkers have pro- 
posed that a simple ID-elastic foundation model can be 
used to accurately describe the contact between elas- 
tic solids [7H9(. In a recent publication they calculated 
the normal stiffness between elastic solids with ran- 
domly rough but nominally flat surfaces, and argued that 
the results are in good agreement with exact numerical 




FIG. 1: In the elastic foundation model the elastic solid is 
replaced by an array of independent springs. 



results [9]. In this note we will show that in fact the model 
of Popov et al fails even qualitatively to describe the con- 
tact mechanics correctly. 

2 Area of real contact and the average interfa- 
cial separation 

Consider two elastic solids with nominally flat surfaces 
squeezed together by the nominal pressure p = F/Aq. 
The average interfacial separation is denoted by u. As 
p increases, the average interfacial separation u mono- 
tonically decreases, while the area of real contact A in- 
creases U 13 • In earlier publications 11, 14 - 17 it has 
been shown that in a large pressure range A ~ p and 
u ~ Inp. This can be understood as follows: As the load 
increases, existing contact patches grow and new, small 
contacts are formed. This happens in such a way that the 
distribution of contact sizes and local pressures remains 
approximately constant over a wide range of loads [111 Il5| . 
It follows that A ~ p and that the elastic deformation en- 
ergy (per unit nominal contact area), U e \, stored at the 
interface must be proportional to the load or the nominal 
contact pressure [14j: 



U c \ = u p, 



(1) 



where uq is a length parameter of order the root-mean- 
square roughness h rms . Since the elastic energy is equal 
to the work done by the external load (assuming hard- 
wall interactions and no adhesion), it follows that 



P = — 



du 



Combining this with (1) gives 

p = p exp(-u/M ) 



(2) 



where po is an integration constant. The theory of Pers- 
son predicts that uq = ah lms and po = j3E*, where E* is 
the effective elastic modulus and a (of order unity) and 
/3 are dimensionless. Both a and _B on ly depe nd on the 
spectral properties of the surfaced, [3, EB43|- 



In the same pressure range where (2) is valid, the area 
of real contact 
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K p 



(3) 



where £ = (( Vh) 2 ) 1 ! 2 is the surface rms-slope and kb2. 

Eqs. (2) and (3) are only valid at such high pressures 
that multi-asperity contact occurs. At very low pres- 
sures the solids will only make contact in the vicinity of 
the highest asperity. In this finite-size pressure region 
the relation between u and p will exhibit large fluctua- 
tions from one surface realization to another [20(. In the 
study presented below the finite-size region is too small 
to be observed on the linear pressure scale used in Fig. 
2J In Ref. [21, 22 1 we have studied numerically and ana- 
lytically the relation between the interfacial stiffness and 
the squeezing pressure in both the finite size pressure re- 
gion and for higher pressures, and in Sec. 5 we compared 
the results for the stiffness with the ID-elastic foundation 
model of Popov et al. 
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FIG. 2: The surface roughness power spectrum as a function 
of the wavevector (log-log scale) for a self-affine fractal surface 
with a roll-off. 



3 Numerical results for A(p) and u(p) 
In the reduction of dimensionality approach a 3D- 
contact problem is mapped on a ID-elastic foundation 
problem. Here we are interested in the contact between 
two nominally flat but randomly rough surfaces. For fric- 
tionless contact, this problem can be mapped on an elas- 
tic half space with a randomly rough surface in contact 
with a rigid substrate with a flat surface. In the con- 
tact mechanics theory of Popov et al. the roughness of 
the ID-substrate has a power spectrum related to that 
of the original via the equation: 



Cm(q) = nqC 2 T>(q). 



(4) 



The spring constant of the elastic foundation is related to 
the effective (or combined) elastic modulus via k = aE* , 
where a is the spacing between the springs. Using stan- 
dard procedures we have generated randomly rough ID- 
surfaces with the power spectra given by (4). As in an 



earlier study |22J, the original 2D surface is self afline frac- 
tal with the Hurst exponent H = 0.7 (or fractal dimension 
Df = 3 - H = 2.3) and with the small and large cut-off 
wavevectors <?o = 1 an d <7i = 8192. We consider two cases, 
namely when the substrate surface is fractal-like in the 
whole interval go < q < Qi , and when there is a roll-off at 
q r = 8, see Fig. [2 

The red and blue solid lines in Fig. 3 have been calcu- 
lated following the procedure outlined by Popov et al.[23- 
[26l | : Each independent spring is compressed into compli- 
ance by the ID rough surface profile where the profile 
overlaps with the initial relaxed spring positions. In each 
step we calculated the force F and area of contact A. 
The applied force F was calculated as the sum of forces 
for all springs in contact: 



F = kY i Aa 



(5) 



where n is the number of springs in contact, Axj is the 
spring compression^ After this we have calculated the 



area of contact A 24-2 



A = ~Y,(am) 2 , 

i=l 



(6) 



where n c is the number of connected regions. In this 
case all springs in connected regions must be in contact, 
rii is the number of springs in each region, arii are the 
diameters of these regions. In the case of full contact the 
area of contact A = Aq, where 



A 



^(aN) 2 , 



(7) 



where iV is the full number of springs. Then there is only 
a single connected region, the diameter of which is equal 
to the length of the system aN. Using §5§ and ([7]) we 
can also calculate the squeezing pressure p = F/A . The 
interfacial separation u was calculated using the formula: 



1 N 



N 



(8) 



where itj is the distance between the end of each spring 
and the substrate surface. For springs in contact itj = 0. 

Using this simple procedure we have calculated the 
dependencies shown in Fig. 3 by the solid lines. This 
dependencies were calculated in range of coordinates of 
elastic foundation (array of springs) from first contact to 
full contact. 

We will compare the predictions of the theory of Popov 
et al. with numerical exact results obtained as described 
in Ref. 



27 



28j . In brief, this method computes the sur- 
face displacements using a Fourier-transform technique 
with a linear surface Green's function that corresponds 
to Poisson ratio v = 1/2. The interaction with the rigid 
surface is treated as a hard-wall repulsion. 
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FIG. 3: (Color online) (a) The area of real contact A in units 
of the nominal contact area Aq, and (b) the average interfacial 
separation u in units of the rms roughness h lma , as a function 
of the squeezing pressure p in units of the effective elastic 
modulus E* . For self-affine fractal surfaces with H = 0.7 
and rms-slope 0.1. The surfaces have the small and large 
wavevector cut-off qo = 1 and qi = 8192, respectively, and the 
roll-off wavevector q r = 1 (blue curves) and q T = 8 (red curves). 



In Fig. I^a,) we show the calculated normalized contact 
area A/Aq as a function of the squeezing pressure. The 
red and blue squares are the result of a numerical exact 
study and the dashed line the prediction using the theory 
of Persson. Since the two surfaces have the same rms- 
slope the theory predict the same curve for both cases 
which agree well with the numerically exact results. The 
red and blue solid lines are the predictions of the theory 
of Popov et al. Since A(p) approach A much faster 
in the model by Popov et al. then in the numerically 
exact theory, the interfacial stiffness K = -dp/du will 
approach infinity much faster (with increasing pressure) 
in the theory of Popov et al., as compared to our exact 
numerical study. Thus the stiffness relation K(p) will 
also be incorrectly given by the theory of Popov et al. 
(see also Sec. 5). 

Fig. E^b) shows the logarithm of the average interfa- 
cial separation u as a function of the squeezing pressure 
p. Again there is good agreement between the numeri- 
cally exact results and the theory of Persson, while the 
approach of Popov et al. fail qualitatively. The results 
presented using the theory of Popov et al. are obtained 



by averaging the calculated quantities over 100 realiza- 
tions of the rough-line topography, with the ID-power 
spectrum given by (4). In each realization the elastic 
foundation has 8192 springs. 

4 Contact stiffness 

Consider two elastic solids with nominally flat surfaces 
squeezed together by the nominal pressure p = F/A . 
From @ it follows that the contact stiffness 

K - d P- P 
du ah lnls 



or 



KK 



E* 



OLE* 



(9) 



This equation is only valid at such high pressures that 
multi-asperity contact occurs. At very low pressures the 
solids will only make contact in the vicinity of the highest 
asperity. In this finite-size pressure region the relation 
between K and p will exhibit large fluctuations from one 
surface realization to another [201] . 

In Ref. [2l| two of us has derived an (approximate) 
analytical expression for the (ensemble averaged) inter- 
facial stiffness in the finite-size region. The derivation as- 
sumes a self-affine fractal surface with the surface rough- 
ness power spectrum shown in Fig. [5J The surface is 
characterized by the Hurst exponent H and the small and 
large wavevector cut-off qo and q\ , as well as a roll-off q r 
(see Fig. [2]). For this model the stiffness per unit area 
in the low pressure, finite size, region is approximately 
given by (see Ref. [2l|) 



E * hki + h) p 
\ L 2 q r I \ h lms I 



K 



1/(1+H) 



■ p 



1/(1+H) 



where L « 27r/go is the linear size of the studied system. 
Note that we can also write this equation as 



Kh r 



\L 2 qJ \E*) 



H/(l+H) 



E* 



(10) 



5 Numerical results for K(p) 

We consider again two cases, when the substrate sur- 
face is fractal- like in the whole interval qo < q < qi, and 
when there is a roll-off at q r = 8, see Fig. [2] In Fig. 
|U[a) we show the calculated interfacial stiffness using the 
theory of Persson. We have plotted Kh Ilns /E* as a func- 
tion of p/E* as these dimensionless quantities enter in 
the theory [see (9) and (10)]. The results in Fig. g^a) 
are in excellent agreement with exact numerical simula- 
tions for the same system, see Fig. 1 in Ref. 22]. One 
can distinguish three regions in the stiffness K(p) rela- 
tion. For very small pressures the stiffness increases as 
K ~ p 1 /( 1+H \ This is a finite size effect, which occurs 
when a single effective Hertz contact region, formed at 
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FIG. 4: (Color online) Log-log plot of the nondimensional 
contact stiffness Kh Ilns /E* vs. nondimensional pressure p/E* 
for self-affine fractal surfaces with H = 0.7 and rms slope 0.1. 
The surfaces have the small and large wavevector cut-off qo = 1 
and qi = 8192, respectively, and the roll-off wavevector q r = 1 
(blue curves) and q r = 8 (red curves). The result in (a) is from 
the Persson contact mechanics theory, which agrees with the 
exact numerical study presented in Ref. The result in 

(b) is using the theory of Popov et al. The vertical dashed 
lines indicate the pressures where the Popov et al. theory 
starts to deviate from the analytical results. 



the highest substrate asperity, prevails. For higher pres- 
sures a region where K ~ p is observed. This region, 
which becomes wider as the width of the roll-off region 
increases, results from contact with many asperities, and 
depends crucially on the long-range elastic coupling be- 
tween the contact regions. Finally, for very large pressure 
the interfacial separation approaches zero and the inter- 
facial stiffness increases towards infinite. 

Fig. QJb) shows the results using the contact mechanics 
theory of Popov et al. The results are obtained by av- 
eraging the contact stiffness obtained in 100 realizations 
of the rough-line topography, with the ID-power spec- 
trum given by (4). In each realization the elastic foun- 
dation has 8192 springs. Fig. |4](b) shows that the the- 
ory correctly predicts the initial (low pressure) relation 
K ~ This result is expected because the study 



in Ref. |2lj shows that K ~ p 1 /( 1+H '> holds even when one 
neglect the elastic coupling between the asperity contact 
regions. However, the region where K increases linear 



with the pressure p is absent in Fig. @|b). This is also 
expected because the K ~ p result depends crucially on 
the elastic coupling between the asperity contact regions, 
which is not included in the theory of Popov et al. As 
shown in Fig. |U[a), the linear region is particular large 
when there is a roll-off in the power spectrum. Most sur- 
faces of engineering interest exhibit a roll-off even larger 
than for the q T /qo = 8 case shown in Fig. [2] Thus in most 
practical applications, in particular involving elastically 
soft materials like rubber, one will be in the linear K ~ p 
region where the approach of Popov et al fails. 

We note that wether there is a roll-off or a cut-off at 
q T has very little influence on the result. However, this 
does not imply that the only thing which matters is the 
range over which the surface is self-affine fractal. The 
point is that including a roll-off or cut-off at q r > qo imp- 
ies roughly that the surface is "periodically" repeated 
(<7r/<7o) 2 times. This implies that there will be many as- 
perities of similar height as the highest asperity. This 
in turn means that the contact will much quicker (with 
increasing pressure) come into the multi-asperity contact 
configuration where the stiffness K depends linearly on 
the nominal squeezing pressure p. This is the origin of 
why the linear relation between K and p starts at lower 
pressures, and extends over a larger pressure range, when 
the surface has a roll-off or cut-off. 

The vertical dashed lines in Fig. 0] indicate the pres- 
sures where the theory of Popov et al. starts to deviate 
from the analytical results in Fig. Ufa). Note that these 
points correspond to the start of the linear K ~ p region 
in the analytical theory. This is expected as the lin- 
ear region corresponds to multi- asperity contact, where 
the elastic coupling between the asperity contact regions, 
which is neglected in the Popov et al. theory, becomes 
important. Thus it is the neglect of the long range elastic- 
ity, and not the reduction in dimensionality, which is the 
basic problem with the approach of Popov et al. In fact, 
a recent study by Scaraggi et al. [29| has shown that if the 
long-range elastic coupling is included in the analysis, it 
is possible to make 2D isotropic roughness approximately 
equivalent to ID roughness. 

6 Summary and conclusion 

We have presented a detailed comparison of the theory 
of Popov et al. with numerical exact results and analyti- 
cal results for self-affine fractal surfaces with and without 
a roll-off. The theory of Popov et al. fails qualitatively to 
describe the A(p) and u(p) relations, and we attribute 
this to the absence of the elastic coupling between the 
asperity contact regions in the approach of Popov et al. 
In fact, a recent study by Scaraggi et ahfiij has shown 
that if the long-range elastic coupling is included in the 
analysis, it is possible to make 2D isotropic roughness 
approximately equivalent to ID roughness. 

We have also presented a detailed comparison of the 
theory of Popov et al. with analytical results for the nor- 



mal contact stiffness. For the case of a roll-off at q r = 8 the 
Persson theory and the exact numerical results presented 
in Ref . [22I ] , exhibit a linear K ~ p region extending over 
3 decades in pressure, while there is no linear region in 
the theory of Popov ct al. The latter theory predicts 
K ~ in the limit of small pressures, but this re- 

sult is expected since an effective Hertz single-asperity 
contact prevails in this case (see Ref. (HJ). However, 
since the elastic coupling between the asperity contact 
regions is absent in the approach of Popov et al., no lin- 
ear K ~ p region is obtained. 
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